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Abstract 


The partitioning of energy and the distribution of the resultant 
ejecta on the moon is numerically modeled using a Eulerian finite- 
difference grid. The impact of an iron meteoroid at 15 km/sec on a gabbroic 
anorthosite lunar crust is examined in detail. The high-speed impact- 
induced flow is described over the entire hydrodynamic regime from a time 
where the peak pressures are 6 Mbar until the stresses everywhere in the 
flow are linearly elastic, and less than 5 kbar. For 5-cm radius projectile 
the latter condition is achieved some o*0.5 msec after impact. The effect 
of taking into account the shock- induced polymorphic phase changes, in 
the plagioclase and pyroxene structure (in gabbro) to the hollandite and 
perovskite structures, respectively, and the subsequent reversion to 
low-pressure phases is demonstrated to enhance shock-wave attenuation. 

A rate-dependent equation of state, is used for describing the hysteretic 
effect of the phase change. The ballistic equations for a spherical planet 

taking into account the decrease of gravity with height, are systematically 
applied to material with net velocity away from the moon. The mass of material 
escaping the moon corresponds to some 28% of the mass of meteorite, less than 
previous estimates, and most of the material lost, is lunar crust. Only 0.2% 
of the meteoroid escapes the moon, all in the vapor phase. In the case of 
accreting planets, a relatively sharp decrease in energy and mass lost from 
such impacts occurs when the escape velocities begin to exceed 1 kn/sec. 

This implies that since the fraction of kinetic energy lost is less than 5%, 
Impact heating of lunar-sized planets in the latter stages of accretion is 
efficient. Most of the impact energy remains on the lunar crust (86.1%) 
however the bulk of the impact energy (66.8%) resides in crustal impact 


ejecta which is distributed with a surface density (mass/area) which decays 
-2 715 

as R , where R, is the radius from the impact. At large distances 

2 

from the impact, the ratio of lunar to meteorite ejecta is $»10 , implying 
that the higher concentrations of meteorite components observed in the 
Apollo 17 breccias, resulted from mostly local impacts. 
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INTRODUCTION 

A central problem in describing the evolution of the terrestrial planets, 
especially in their later stages of accretion, is determining the partitioning 
of ejecta, in terms of mass and energy, upon impact of a meteorite. Although 
it is not yet clear whether the terrestrial planets formed from small parti- 
cles in chemical and thermal equilibrium with the primordial gaseous nebulae 
as proposed by Lewis (1972), or from a chemically heterogeneous meteorite- 
like population of larger objects (Urey, 1952) or possibly, inhomogeneously 
directly from a condensing gas cloud (Clark et al. , 1972), the present paper 

is concerned with how a lunar-sized object responds to the impact of meteoroids 

2 

of 1(T jv km-radius. The present surface bears strong evidence that such 
meteoroids have impacted the moon since it evolved to approximately its 
present radius. Since some 80% of the highland samples returned from the 
lunar highlands are now known to be impact breccias with various stress and 
thermal histories, the thermal state and provenance of the ejecta from a 
single hypervelocity impact, needs to be studied in detail. 

Recently, Overbeck et al. (1975) have argued that the ejecta from large 
impacts on the moon give rise to major transport and mixing processes which 
operated on both previous fragmental materials deposited on the moon (the 
regolith), as well as presumably the initial ejecta from the primary impact 
itself. The available evidence for such processes which may play an important 
role in modifying cratered surfaces, are based on photogeology and terrestrial 
cratering studies and subsequent theoretical analysis. This phenomena is not 
described in the present work, however. 

Previously (O'Keefe and Ahrens, 1975, 1976) have examined in some detail 
the question of partitioning of energy from a large impact on the moon and 
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have, with assumptions regarding the properties of the meteorite and the 
lunar crust, calculated the decay, in space and time, of the strong shock 
which is induced by the impact. The shock is expected to alter the minerals 
of the crust as well as provide melted rock, unlike that resulting from 
normal igneous processes (Warner et al . , 1973; Simonds, 1975). 

Previous knowledge of ejecta distributions from hypervelocity impacts 
onto geologic materials is confined to photographic recording of the 
formation of the ejecta plume, transient cavity, and post-impact analyses of 
mass movement. Impact flows from mm-sized projectiles at speeds of 6-7 km/sec 
were studied. These experiments employed quartz sand (Braslau, 1970 and 
Stbffler et al. , 1975) and solid basalt (Gault et al. , 1963) target media 
and have provided considerable data on the general impact energy partitioning, 
and the shape and evolution of the ejecta plume. Although carried out under 
the earth's gravity field at pressures of 3 to 0.5 mm Hg, the quantitative 
post-impact ejecta distributions measured by Stbffler et al . , have provided 
some experimental support for the relatively uncertain ejecta thickness 
distributions assumed by McGetchin et al. (1973) to infer thicknesses of basin 
deposits at the various Apollo landing sites. Although the quantity of 
meteoroid material in the present calculation is small, (-'>4 kg), and since 
we have not included gravity body forces in the cratering flow itself, we 
believe the results obtained in the large Froude number regime can be scaled 
to large Impacts to the degree that the moon's curvature in the region of 
high-pressure flow can be neglected. 

The objective of the present study is to examine the partitioning of 
impact energy, into the kinetic and internal energy of the ejecta, and to 
examine the spatial distribution of ejecta on the moon. The spatial 
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distribution of the meteorite debris and lunar-surface ejecta are separately 
treated. We also consider briefly the implications of this, initial, calcu- 
lation on the larger question of accretionary processes for the terrestrial 
planets. 

CALCULATIONAL OBJECTIVES 

The present results were obtained by using the conservation of mass, 
momentum, and energy equations in finite-difference form (Hageman and Walsh, 
1970) within a Eulerian framework with approximate equations-of-state as 
outlined previously (O'Keefe and Ahrens, 1975). 

For this initial calculation it was assumed, that the meteorite had an 
iron-like composition because the appropriate equation-of-state is well 
known. Although differentiating the meteorite debris from lunar surface 
fragments in the calculation turns out not to be a problem, initially it 
presented conceptual difficulties. Unfortunately, an iron meteoroid in 
some respects, is probably atypical of objects which have impacted the moon, 
in light of the gross depletion of iron on, or within, the moon (compared to 
the other terrestrial planets) . With this drawback in mind in the present 
calculation, the impact of an iron object is assumed at 15 km/sec. This 
velocity lies within an average range of from 15 to 18 km/sec inferred by 
Zook (1975) for the present moon and is similar to meteoroid velocities 
observed as radiant shows by the Prairie Network (Wassen, 1973, p. 137). 
Impact velocities on the order of at least ~ 10 km/sec can. Independently, 
also be inferred on the basis of the present knowledge of the Hugoniot curves 
for rocks, and the shock-metamorphlc effects, especially melting, associated 
with rocks sudied in-3itu and within the ejecta, from many of the impact 
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craters found on the earth. The present calculation therefore deals with 
the ejecta from the flow produced by a 15 km/ sec Impact of an Iron object 
on a gabbrolc anorthosite half-space (In the vicinity of the impact, ~ 100 
meteorite radii). However, once the initial ejecta velocity spectra are 
calculated, the ejecta paths are calculated for a spherical moon. The 
equation of state of the moon is modeled using data from sample 15,418 
(Ahrens et al. , 1973). 

Two new features now incorporated in the calculation which should make 
the present work more credible and useful are: 

(1) An explicit description of the kinetics of shock-induced phase 
changes for silicates comprising the lunar crust. This is provided by 
an algorithm which in essence prescribes a time-dependent equation-of-state. 
This brings the description of silicates into agreement with the release 
adiabat data, obtained in the mixed-phase regime reported by Ahrens et al. 
(1969), Petersen et_ al. (1970) and Grady et al. (1974). 

(2) Inclusion of the necessary ballistic equations to calculate the 
trajectories of ejecta. Of specific interest are the mass and internal and 
kinetic energy content of both the initial ejecta which escapes the moon, 
or another terrestrial planet, as a function of escape velocity, as well as 
the spatial distribution of ejecta density at various times on a spherical 
moon. 


EQUATIONS OF STATE 

The shock wave data for plagioclase and pyroxene (McQueen at al . , 1967; 
Ahrens et al. , 1969; and Ahrens and Gaffney, 1972) demonstrate that both the 
component minerals, as well as rocks containing these minerals, undergo a 
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probable succession of shock-induced changes above ~ 125 to 150 kbars. The 
occurrence of these solid-solid transitions is not a trivial detail, since 
they account for ~ 351 increase in equivalent zero-pressure density. This 
results in almost an order of magnitude more compression, than the intrinsic 
self-compression of the minerals up to the transition pressure! Complete 
transformation to high-pressure phases is expected to occur in the 500 to 
700 kbar range in a lunar crustal rock such as 15,418 (Ahrens et al. , 1973)* 
Coincidentally upon release of shock pressure from this pressure range the 
onset of melting occurs (De Carli and Milton, 1965 and Kief fer et al. , 1976) . 
Both the Huge dot data and the results of static high-pressure phase studies 
demonstrate that the high-pressure phases most probably produced in the shock 
case for plagiodase is in the hollandite structure (Ahrens et al. , 1969) and 
for pyroxene in the perovskite structure (Liu, 1975). The latter result was 
anticipated in interpreting the initial results for gabbroic anorthosite 
(Ahrens et al . , 1973). Taking this into account, we have separately fitted 
an equation-of-state of the Tillotson form for each phase (O'Keefe and Ahrens, 
1975). In the compressed region (p > p q ) the pressure is given by 

p - | a + b/[E/(E o 

where P is pressure, p q is the initial density, E specific internal energy, 

n = p/p and y » n- 1. Here (a + b) and. A, are the zero-pressure and low- 
o 

temperature, Gruneisen parameter and bulk modulus, respectively. The 
constant, a, was taken to be 0.5 so that at high-energy densities the 
thermal pressure would approach that of an electron gas. For low densities 
and high temperatures, such that p < P q and with E > E^, the form used is: 


n 2 ) 


+ 1] J Ep + Ay + By 


( 1 ) 
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where a end 8 were chosen so Chet the equation approaches a polytropic 
equatlon-of-state at low densities and high specific internal energies with 
a polytropic exponent of (a + 1). 

In Table 1, the parameters for the high- and low-pressure phases 
(Fig. 1) are given, and E and E' represent the energy for incipient and 
complete vaporization at atmospheric pressure. In order to achieve unloading 
along isentropes which approximate those observed for feldspars by Ahrens 
et al. (1969) and Grady et al . (1975) above ~150 kbar, we have employed the 
phenomenological model of Horie (1966). The total energy and specific volume 
are given by 


f V, (P.E) 
eq hpp 


(1-f ) 

eq 


V. (P,E) 
ipp 


(3) 


E 


f E, (P.E) +(l-r ) E. (V,E) 

eq hpp eq Jtpp 


(4) 


where the subscripts hpp and ipp indicate low- and high-pressure phases, and 
thermodynamic equilibrium in the mixed phase region is treated the same as 
in O'Keefe and Ahrens (1976). The time variation of the mass fraction, f, 
of material in the high-pressure phase is assumed to be given by 


df/dt - (f -f)/t (5) 

eq 

where f is the equilibrium (minimum free-energy) mass fraction and t is 
eq 

the time constant of reaction. As indicated in Table 1, the increase in 
internal energy upon transforming from the low- to high-pressure phase is 
taken as 0.013 Mbar-cm /g, independent of temperature. Based on the 
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experimental unloading profiles of Grady et al. (1974) for release data 

taken in the quartz-stlshovlte regime, t, is assumed to be ~ 0.25 Msec. 

It is worth commenting on the consistency of the constants of Table 1, 

based largely on Hugonlot data and Thomas- Fermi calculations, with other 

more conventional shock wave reductions of data, and the high-temperature, 

low-pressure properties. This can sensibly be done only in the case of iron, 

for which the present results predict incipient and complete melting to occur 

upon ad. stbatic release from principal Hugoniot states at ~ 0.63 and ~ 0.93 

Mbar, respectively. Whereas, post-shock temperatures comparable to the melting 

point (1808 # K) are calculated, to occur upon pressure release from ~ 1.7 Mbar 

(McQueen et si., 1970). Moreover the Tiilotson formulation (Eq. 1.) over- 

3 

estimates the post-shock, complete-melt densifv at 7.39 g/cm , whereas, the 

3 

measured density of liquid iron, at the melting point is <■« 7.00 to 7.01 g/cm . 
Finally, conventional Mie-Gruneisen reduction of the shock wave data predicts 

3 

a post-shock density, which is too low, 6.704 g/cm , at the melting point 
(McQueen <it al . , 1970). We suspect some of these discrepancies may arise 
from a lack of properly taking into account the a-e phase change and what is 
now known about the iron melting relations (Liu and Bassett, 1975). 

Since the present high-speed flow calculations can be carried out to 
relatively long times at which point the stress field becomes linearly 
elastic, the yield criteria and complete elastic constants arising from an 
assumed Poisson's ratio of 0.20 must be utilized. 

The finite-difference formulation permits a full tensorial treatment in 
cylindrical coordinates of stress and strain in the finite strength regime. 

A simple model which will account for dynamic yielding should, for example, 
specify the maximum shear-stress "or strength" supportable by the medium. 
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allow this strength to Increase with Increasing mean principal stress, and 

also account for the thermal weakening of the material by allowing the 

yield strength, Y, to decrease with increasing internal energy content at 

2 

constant strain. Using the second stress invariant, 2Y , as a measure of 
the yield strength, following the treament of Hageman and Walsh (1970), we 
assumed the following form for Y for gabbroic anorthosite: 

Y(kbar) - (2.69 + 337. 6y - 901 y 2 ) (1 - E/E ) (6) 

m 

where E is the internal energy density required for incipient melting under 

IB 

standard conditions. Thus, in the calculation when the second stress 

2 

invariant exceeds ?Y all stresses are reduced proportionately so as not 
to exceed the yield strength specified by Eq. 6. 

The above Mohr-Coulomb type yield condition is based on at present, a 
very limited set of data for the dynamic yielding of plagioclase- and 
pyroxene-bearing rock. At zero-confining stress. Eq. 6 satisfies the failure 
data of Kumar (1968) for one-dimensional stress loading of a terrestrial 
basalt at rates of ~ 10^ kbar/sec and the Hugoniot-slastlc-limlt data of 
Ahrens et al. (1973) at loading rates of ~ lO 1 ^ kbar/sec for lunar gabbroic 
anorthosite (sample 15,418). Finally, although the Mohr-Coulomb criterion 
is applied for dynamic yielding under compressive loading, we have utilized 
the value of 0.5 kbar for the dynamic tensile strength of gaboroic anorthosite. 
(This value is certain to be dependent on rock volatile content.) Although 
some data on dynamic tensile failure for rocks are available (Shockley et al. , 
1975) we suspect our description of the latter stages of crater excavation, 
and hence, ejecta distribution within several crater radii, are sensitive to 
this poorly constrained quantity. 
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THE EFFECT OF THE SHOCK- INDUCED Piv..E CHANGE 
ON THE IMPACT FLOW 

Sine* previous calculations of impacts iu silicate material have used 
equation-of- state formulations which simply faired through both the high- 
and low-pressure portions of the Hugoniot (O'Keefe and Ahrens, 1975; Bjork, 
1956; and Gault and Heitowit, 1963), it i3 interesting to examine the 
details of the early crateiing flow which ere affected by the shock- induced 
phase change in the silicate. 

In the present calculations, a spherical iron meteoroid was assumed to 
impact a gabbroic anorthosita lunar surface at a velocity of 15 km/ sec, this 
calculation is similar in many respects to that previously described (O'Keefe 
and Ahrens, 1975). However, previously the radius of the meteoroid was 23.2 
km and in this calculation the radius is 5 cm. The size of impacting body 
was made relatively small, because we wish to include the phase change 
kinetics and also because of considerations of numerical viscosity. To 
obtain the necessary numerical resolution so as to examine phase change 
evolution in time and space, the time step which is a function of the sonic 
plus mass propagation time across a cell, must be much less than the 
characteristic phase transformation time (t ■ 0.25 usee). This criterion 
restricts the size of the impacting object for acceptable cell numbers 
(<: 3000) to radii less than ~10 cm. In addition, the numerical viscosity 

is less for smaller mesh spacing. 

We want to emphasize that even though the present calculation is for a 
5 cm-radius impacting body, these results can be scaled. This impact calcu- 
lation is valid at early times (when the mean particle velocity is greater 
than the escape velocity from the crater being excavated) for all bodies 


12 . 


whose radii (r) are ouch greater than the sonic velocity (c) times the 
characteristic phase-change time (r >> tc); in other words, these bodies 
for which the shock propagation time across the body Is much greater than 
the characteristic phase-change time. The flow field obtained In the present 
calculation can be scaled to that for a nominal Mare Imbrium-type Impact 
described by O'Keefe and Ahrens (1975) in which a 23.2 km-radlus iron 
projectile is assumed as the impactor. The distance and time are scaled 
(hydrodynamically) using the ratio of projectile radii, r^/r^ “ 4.64 x 10^. 

Thus for a time, t 2 ■ 25.4 psec, for the present problem, the equivalent 
time in th~ nitial report of (O'Keefe and Anrens, 1975) is ~ 11.8 sec. 

The flow fields associated with both calculations, having a similar peak 
pressure, are compared in Fig. 2 and 3 and some key parameters of the two 
flows are listed in Table 2. In both flow fields the peak pressure is 
~ 0.66 Mbar, however, the maximum depth of penetration of the meteorite 
without a phase change io 65 km compared to 54 km with a phase change. The 
most significant difference between the two flows is the rapid attenuation 
of shock pressure when the phase change is taken into account. Thin is 
demonstrated in Table 2. In our earlier calculation, the peak pressure 
extends to a depth of Ilf km versus 66 km for the present case. The rapid 
attenuation results from the high rarefaction-wave velocity associated with 
the high-pressure phase material. Earlier, at this stage in the flow, (Fig. 2) 
the meteorite lay completely below the initial lunar surface, whereas in 
Fig. 3 a portion of the vaporized iron object is 5 km above the initial 
surface. Moreover, comparison of Figs. 2d and 3d demonstrates that much of 
this vaporized meteorite material is moving inward when phase changes are 
taken into account, whereas before the entire flow is still radially outward 
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at this time. 


EJECTA BALLISTICS 


At a series of times, in the present calculation, we have employed the 
ballistic equations for describing ejecta for which the net particle velocity 
of the flow is pointed above the local, flat, horizon. Taking the sphericity 
of the moon, and the decrease in gravity with height into account, yields 
for the range (Thomson, 1963 


R 2R tan 
o 


-1 


—2 

V sin 8 cos g 
o 

(1-V 2 cos 2 6) 


(7) 


where 



o 


2V 2 /R g 
o oo 


( 8 ) 


Here R q is the mean lunar radius, 1738 km, V q is the magnitude of the 

velocity vector of the contents of a eulerian cell , g the surface 

2 

gravitational acceleration, 161.8 cm/ sec , and 8 is the elevation angle of 
the particle, from the local horizontal. 

EJECTA AND PLANETARY ACCRETION 

The relative amount of mass exceeding various escape velocities as a 
function of time is shown in Figure 4 for times early in the flow. It is 
significant that the amount of mass exceeding the escape velocity of the moon 
is approximately 20 percent of the mass of meteoroid. 

At a much later time, 545 ysec, the amount of rock escaping the moon 
has increased only slightly, ~ 28% of the mass of the meteoroid, however. 
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only ~ 2% of the escaping lunar material is vaporized as a result of the 
impact. Host of the meteoroid remains on the moon, our present calculations 
indicate, only 0.2% escapes, all in the vapor phase. 

The present amount ot material escaping the moon is less than would be 
predicted from the results of Gault ejc al . (1953) which are based on experi- 
ments with aluminum projectiles having diameters less t ■ an 3 mm. Moreover, 
they infer on the basis of high-speed framing camera photography that for 
impacts of 6.1 to 6.4 km/ sec on the order of an equivalent mass of the 
projectile, probably largely target derived material, is ejected at speeds 
greater than the lunar escape velocity. Note that in the present calculations 
the mass exce* !,ng the escape velocity of the moon reached its final value 
early in the impact flow (within 5 ysec, or in terms of larger-scale events, 
a time interval corresponding to translation of the meteoroid through a 
distance about equal to its radius). However, the amount of mass exceeding 
lower escape velocities, e.g., velocities less than ~ 1 km/sec, is still 
increasing at the end of/ysec. Referring to Fig. 5, the material exceeding 
an escape velocity of 1 km/sec is primarily in the liquid and vapor state, 
whereas the material having escape velocities less than 1 km/sec is primarily 
in the solid state, and the amount in the solid state is still increasing at 
the end of 25 ysec. One interesting feature of this figure is the sharp 
cutoff in mass lost at escape velocities slightly exceeding 1 km/sec. The 
implication is that for 15 km/sec meteoroids, the moon and larger bodies are 
efficient accretors, whereas for smaller bodies having escape velocities much 
less than a kilometer per second, the amount of mass lost exceeds the mass 
of the impacting meteoroid. The relative amount of energy lost as a function 
of escape velocity is shown in ire 6. For impacts on bodies having 
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escape velocities greater than 1 km/ sec, the fraction of energy lost is 
less than 5% and therefore impact heating of terrestrial planets in the 
latter stages of accretion is an efficient process. 

Finally, we have examined the ejecta distribution at late times 
( ~ 545 ysec), both to provide a first-order determination of crater volume, 
and secondly, to examine the ejecta blanket height versus radius in light 
of both the small-scale laboratory results, and the explosive data which 
have already been applied to the moon by McGetchin et al. (1972). The 
mass-flux distribution as a function of angle of ejection and velocity 
weighted by mass versus angle of ejection are shown in Figs. 7 and 8. This 
calculation is again for a 5-cm radius iron projectile (4 kg), impacting the 
moon at 15 km/sec. The shock pressures present in the impact- induced flow 
range from nearly 2 Mbar, at 9 ysec to ~ 5 kbar, at 545 ysec. The ejecta 
plume is undoubtedly rather smooth versus the irregular pattern calculated, 
which demonstrates the effects of finite-sized zones. Even with 167 cells 
used to describe the meteorite at initial impact, the zoning used is too 
coarse to resolve the jetting which is known to occur under oblique conditions 
(Kieffer, 1975). The mass-weighted velocity versus angle plotted in Fig. 8, 
is effectively a plot of momentum distribution, in that the quantity summed 
per each 10° increment of angle is 

E®! (9 > 

i 

where V q has the same meaning as in Eq. 8 and the summation is carried out 
over all particles with a net axial velocity away from the lunar surface. 

Note that taken together, Figs. 4, 7, and 8, indicate that while high-shock 
pressures are present in the flow, e.g. at 9 ysec, the mass of ejecta 


£ v „i -i 
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present is of the order of the mass of the meteoroid and has a velocity on 
the order of ~ 10"* cm/sec. (This is the only period during the cratering 
process in which material is lost from a moon of present size.) Moreover, 
at successively later times, as the mean pressure decays, so does the mass- 
weighted velocity and the ejecta achieves a mass-weighted velocity which is 
of the order of ^10 to 10 cm/sec. We suspect that the curves shown for 
462 and 545 psec, are dependent on the dynamic tensile criterion used in 
the calculation. 

Comparison of the present ejecta patterns with small-scale laboratory 

experiments, largely carried out at impact exceeds significantly lower than 

those required to produce substantial vapor (Ahrens and O’Keefe, 1972) present 

some interesting contrasts. It should immediately be recognized that the 

expansion of the inner ejecta plume of vaporized meteorite depicted in 

Fig. 3 will give rise to late-time ejecta which have no counterpart in the 

experiments of Gault el: al. 1963; Oberbeck and Morrison, 1976; and Schneider, 

1975. Hence comparison with these results should involve ejecta distributions 

at relatively early times, i.e., 9 and 47 psec, when the interacting shock 

and rarefaction wave effectively launch largely material from near the target 

/ at 

free-surface. Comparison of Figs. 7 and 8 indicate that/early times the 
bulk of the low-speed ejecta is launched at low angles, whereas the higher 
speed ejecta (which is also hottest) is launched at higher angles. Quali- 
tatively this is observed by Overbeck and Morrison in their shock launched 
ejecta. In our Eulerian calculation the entire flow field is in motion 

by the time the stresses are low, we cannot determine the initial position 
of each mass particle in the flow. It is interesting to note that the present 
ejecta flux is far more evenly distributed with respect elevation than would 
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be Inferred from Gault et al. 1963 framing sequences. Like Gault et_ al . 

1963 and McDonnel et al. 1976 (both experiments carried out on crystalline 
silicates) and unlike the results of Schneider 1975 (carried out on a glass), 
the calculated ejecta flux is a maximum (at ~ 47 Msec) at high elevation 
angles. It is also interesting that our calculation qualitatively gives an ejecta 
flux distribution with elevation having a broad minimum at ~ 45° which resembles 
that of McDonnell et al. However, probably because of the gross velocity 
differences and rock strength and surface effects the relative amount of 

3 

ejecta in our calculation is ~ 2 x 10 greater. 

We note that at late times (545 Msec) the total ejecta is about 1.6 x 

g 

10 grams for a meteorite having an initial mass of 4.079 kg. Thus the 
total amount of ejecta is about 4 x 10"* times the mass of the meteorite. 

3 

If a mean density of 2.7g/cm for lunar crust and a crater depth-to-radius 
ratio of 0.4 is assumed for a conical shaped crater, an equivalent radius 
of 521 cm crater is calculated. This is plotted for comparison with ter- 
restrial and lunar craters in Fig. 9. Our theoretical crater is seen to 
be similar in size to those mapped by Moore (1976) and perhaps most comparable 
to the Ranger craters produced by slower, but more massive projectiles 
impacting the more distended regolith. The present calculation gives the 
final mass of ejected meteorite as ~1500 grams, or one third of its initial 
mass. This value is again probably heavily affected by the rheology assumed 
for the meteoroid. A fluid was assumed for the present case. 

Table 3 summarizes the partitioning of energy at late times (545 Msec). 

Of special Interest is the great efficiency of the impact process in trans- 
ferring kinetic energy into the thermal energy at the later stages of 
planetary accretion. As can be seen from the Table most of this thermal 
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energy resides in the lunar-surface ejecta which is spread over the moon in 
accordance with the distribution shown in Fig. 10. The ejecta that travels 
the farthest is the initial ejecta. This is also the hottest. Since the 
stresses at late times have all decayed into the assumed elastic range, the 
present calculation implies that little of the kinetic energy, 1.1% of the 
projectile energy resides in the non-ejecta component. The 1.1% value, 
then represents energy which ultimately will largely remain as internal energy 
in the vicinity of the impact as only ~0.01 to 0.1% of the impact energy is 
inferred to convert to seismic energy (Schultz and Gault, 1975). 

The ejecta distributions for meteoroid fragments and lunar surface 
ejecta for early and late times are shown in Fig. 10. These are obtained 
using the trajectory of the contents of each cell moving away from the lunar 
surface. At 545 Msec there are some 123 cells of variable mass in such 
trajectories, each with a mean internal energy. The ejecta surface density 
is calculated by determining the end point of each trajectory on a spherical 
moon. 

At late times at distances, *>10 cm from the impact, the calculated 
ejecta surface density, A, on a spherical moon can be fit with 

A(g/cm 2 ) « 1.29 x 10^ R 2 ‘ 7 ^ (10) 

where R(cm) has the same meaning as in Eq. 7. Subdividing the ejecta distri- 
bution into close-in, and far-out, regimes is clearly arbitrary, since the 
first-three shortest radii surface densities calculated, lie within the 
hypothetical crater. At late times the meteorite ejecta, which is very 
coarsely zoned (i.e. only three particles!), can be fit with 

A(g/cm 2 ) - 4.42 x 10 7 r “ 2 ' 930 


(ID 
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These observed ejecta densities decay with radii slightly less rapidly 
than the observations made on a half space, at low atmospheric pressures, 

for sand, in the earth's gravity field by Stttffler et al. They observed a 

—3 26 o 

decay of R . McGetchin et al. assumed a value of R * for estimating 

thicknesses of basin deposits on a flat moon. The ratio, T, of lunar rock 

2 

ejecta to meteorite ejecta is >10 . This ratio, to first order, should 
be Independent of projectile size since the ranges calculated are independent 
of projectile mass. The present 

values for T are somewhat greater than the values inferred for Apollo 17 
breccias by Morgan et al. (1974). Taking the Morgan et al data and the 
calculated radius dependence for T at face value implies that if the sidero- 
phile minor element contents measured in the Apollo 17 breccia boulders come 
from tne one projectile, and different breccia samples have geochemical 
imprints of different objects, the impact events occured at ranges where T 
< 10 , or at R £ 1 km, from the Apollo 17 site and hence do not represent 

the major basin forming imp? :ts . 
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Table 1. Equation of State Paraaeters 
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Table 2 


Comparison of 1; vact flow fields for 23.2 ka radius, IS km/s iron meteoroid 
at 9. 35 sac a phase change (Fig. 2) as compared to flow for phase cliange 
in gabhroic anorthosite (Fig. 3) 25. A nsec after impact, distance scaled 
to 23.2 ha-radius impacting object. 


Parameter 

Units 

Fig. 2 

(no phase change) 

Fig. 3 

(phase change) 

Min. depth meteorite 

km 

23 

-3 

Max. depth meteorite 

km 

65 

54 

Max. diameter deformed 
meteorite 

km 

115 

103.5 

Peak pressure 

Mbar 

0.662 

0.66 


Depth from surface to 
11 peak pressure along 

centerline km 56 47 

Depth from surface to 
501 peak pressure 

along centerline km 100 53 

Depth to 1001 peak 
pressure along 

centerline km 118 66 


Peak pressure, meteorite- 

moon interface Mbar 


0.041 0.396 


Max. depth of shock along 
centerline (11 peak 
pressure on leading 
edge) km 


140 


75 



Table 3. Gross Energy Partitioning for 15 km/sec 


Impact of 

Iron Object onto Gabbroic 

Anorthosite 



Kinetic Energy 

Internal Energy 

Total 


(%) 

(%) 

m 

Lunar Surface 

6.9 

86.1 

93.0 

(ejecta)* 

(5.8) 

(66.8) 



Iron Object 

0.0205 

6.7 

6.7 

(ejecta) 

(.0062) 

(2.4) 



* Indicated by parenthesis 



Figure Captions 


Figure 1. 


Figure 2. 


Figure 3. 


Equation of state of gabbroic anorthosite in low (ipp) and 
high pressure (hpp) regimes. Log^g (pressure) versus com- 
pression . The pressure range assumed for mixed-phase regime is 
indicated. Equation of state of 2.86 g/cm basalt used in Gault 
and Heitowit (1963) (G & H) formulation is shown also for 
comparison. In the G & H treatment release adiabats are 
assumed to coincide with Hugoniot curve. 

Hypervelocity flow field for 23.2 km radius, 15 km/sec iron 

meteorite striking a gabbroic anorthosite lunar crust, 9.35 sec 

after impact (scale time for 5-cm radius meteoroid is 20.2 psec). 

Phase changes have not been specifically taken into account. 

Contour values give percentage of difference between minimum 

(min) and maximum (max) value (cgs) indicated for each variable. (a 

9 12 

pressure, min = -2.25 x 10 , max = 6.62 x 10 ; (b) internal 

energy, min = 0.0, max = 1.81 x 10^; (c) axial velocity, min = 

-2.47 x 10^, max = 5."j x 10^; and (d) radial velocity, min = 

0.0, max = 2.81 x lo\ after O’Keefe and Ahrens (1975). Vertical 
and horizontal scales are identical. 

Hypervelocity flow field for 5-cm radius, 15 km/sec iron meteorite 
striking a gabbroic anorthosite lunar crust, 25.4 psec after impact. 
Second scale shown is that appropriate for comparison to 23.2 km 
radius meteorite impact of Fig. 2. Phase changes have been taken 
into account. Contour values give percentage of difference between 
minimum (min) and maximum (max) value (cgs) indicated for each 
variable, (a) pressure, min * -2.0 x 10^, max * 6.6 x 10^; 
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Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Figure 9. 


(b) internal energy, min - -1.3 x 10 10 , max, - 1.2 x 10 11 ; 

(c) axial velocity, min - -2.7 x 10 5 , max * 6.5 x 10 5 ; and 

(d) radial velocity, min ■ -1.9 x 10^, max * 2.4 x 10^. 

, Relative mass escaping planetary gravity field as a function of 
time and escape velocity for 15 km/sec iron meteoroid impacting 
gabbroic anorthosite crust. 

, Lower bound of relative mass partitioned in terms of vapor, liquid, 
and solid escaping planet from 15 km/sec iron meteoroid impact 
as a function of escape velocity at a characteristic time of 
25 ysec for a 5-cm radius meteoroid. 

. Lower bound of relative energy (partioned in terms of that 
contained in vapor, liquid and solid ejecta) escaping planet 
from 15 km/sec iron meteoroid impact as a function of escape 
velocity at a characteristic time of 25 Msec for 5-cm radius 
meteoroid. 

. Ejected mass flux versus ejection angle from lunar surface at 
various times. Flow corresponds to a 5-cm radius projectile 

striking a gabbroic anorthosite target initially at 15 km/ sec. 


Mass weighted, ejecta velocity versus ejection angle from lunar 
surface at various times. (Same flow as in Fig. 7.) 

Logio (crater radius) versus Log^ (total estimated impact 
energy), for various laboratory (Gault, 1973; Fechtig et al. 
1974) and field observations (Moore, 1976, also see O'Keefe 
and Ahrens, 1976). Large symbol indicates theoretical crater 
radius calculated for ejecta mass calculated for 5-cm iron 
projectile impacting gabbroic anorthosite at 15 km/sec. 



Minimum gravitational energy line calculated for bowl-shaped 
3 

crater for, 3 g/cm density rock using formulae of O'Keefe 
and Ahrens, 1976. 

Figure 10. Ejecta density on the moon versus range from center of crater. 

Large symbols are ejecta distribution at 545 usee after impact 
of 15 km/sec, iron meteorite into gabbroic anorthosite. Small 
symbols are distribution of initial, hot, ejecta calculated 
from particle trajectories 25 usee after impact. 
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